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On Quasi-periodic Differential Pencils with Jump Conditions Inside 

the Interval 

V. Yurko 

Abstract. Non-self-adjoint second-order differential pencils on a finite interval with non- 
separated quasi-periodic boundary conditions and jump conditions are studied. We establish 
properties of spectral characteristics and investigate the inverse spectral problem of recover¬ 
ing the operator from its spectral data. For this inverse problem we prove the corresponding 
uniqueness theorem and provide an algorithm for constructing its solution. 

Key words: differential pencils, non-separated boundary conditions, inverse problems 
AMS Classification: 34A55 34B24 34B07 47E05 

1. Introduction 

Consider the boundary value problem B of the form 

y" + (p 2 + PP(x) + q(x))y = 0, x G [0,T], (1) 

y(0) = ay(T), y'(0) - (ipti + h)y(0) = f3y'(T), (2) 

y(bj + 0) = 7 jy{bj - 0), y'(bj + 0) = 7 7 “V(fy - 0) + (i/977' + Vj)y(bj ~ 0), j = 1, N — 1, (3) 

0 = bo < b± < ... < tv_ 1 < b]y = T, 

where p is the spectral parameter, p(x),q(x) are complex-valued functions, and h', h, a, f3, 7, , 
are complex numbers, ^ 0, Zq := a(l=F/r , )+/3 7^ 0, := (7^-f- (7^) —1 )/2=F?7'./2 7^ 

0. Assume that p(x) G AC[0, T] and q(x) E L(0,T). We study a nonlinear inverse problem of 
recovering B from its spectral data. For this inverse problem we prove the uniqueness theorem 
and provide a procedure for constructing its solution. 

Inverse spectral problems often appear in mathematics as well as in applications [1-3]. For 
Sturm-Liouville operators with separated boundary conditions, inverse spectral problems have 
been studied fairly completely (see the monographs [1-3] and the references therein). Such 
problems for Sturm-Liouville operators with non-separated boundary conditions investigated 
in [4-7] and other works. 

Differential equations with nonlinear dependence on the spectral parameter arise in various 
problems of mathematics as well as in applications. In particular, several examples of such 
spectral problems arising in mechanical engineering are provided in the book [8] of Collatz; see 
also [9-11], where further references and links to applications can be found. Detailed studies on 
direct spectral problems for some classes of ordinary differential operators depending nonlinearly 
on the spectral parameter can be found in various publications, see e.g. [9-11]. 

Inverse spectral problems for differential pencils, because of their nonlinearity, are more 
difficult for investigating, and nowadays there are only isolated fragments, not constituting a 
general picture, in the inverse problem theory for equation (1). Some aspects of the inverse 
problem theory for pencils under various restrictions were studied in [12]-[16] and other works, 
but mostly only particular questions are considered there. Inverse problems for general non- 
selfadjoint boundary value problem (l)-(3) with jump conditions inside the interval have not 
been studied yet. We note the inverse problem, considered in this paper, appears in the inverse 
problem theory for differential operators on spatial networks with cycles (see [17-19]) which 
have many applications in natural sciences and engineering. 

Some words about the structure of the paper. The statement of the inverse problem is 
provided in section 2. In Section 3 properties of the spectrum are established. In particular, 
the Weyl-type function and the corresponding Weyl sequence are introduced and investigated. 
In Section 4 we provide the solution of the inverse spectral problem for the boundary value 
problem B. 
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2. Statement of the inverse problem. 


Denote by S(x,p ) and C(x,p ) the solutions of equation (1) satisfying jump conditions 
(3) and the initial conditions 

S( 0, P ) = C'( 0, p) = 0, S'{ 0, p) = C{ 0, p) = 1. 

For each fixed x, the functions S ( ' p \x^p) and C ( ' v \x 1 p ), v — 0,1, are entire in p of expo¬ 
nential type, and (C(x, p), S(x, p)) = 1, where (y,z)\=yz' — y'z is the Wronskian of y and 
2 . Put (p(x, p) = C(x,p) + {iph! + h)S(x, p), d(p) = S(T,p), cfi(p) = C(T,p). Eigenvalues 
V = {p n }ngz of the boundary value problem (l)-(3) coincide with the zeros (counting with 
multiplicities) of the characteristic function 

a(p) = a<p(T, p) + /3S'(T, p) - (1 + a(3). (4) 

Let A := {n : n = ±1,±2,...} = Z \ {0}, and let V = {v n }neA be zeros (counting with 

multiplicities) of d(p). Then {u n } ne a are the eigenvalues of the boundary value problem B 
for Eq. (1) with jump conditions (3) and with the boundary conditions y( 0) = y(T) = 0. The 
function d(p) is called the characteristic function for B . Without loss of generality, we agree 
that the numeration is chosen such that v n ^ if nk < 0. Let m n be the multiplicity of 

u n ( v n = V n + 1 = • • • = Vn +mn -i )• Put /:={n6 A : V n - 1 Vn}, I' := {n G I : m n > 1}. 

Denote D(p) = a<p(T, p) + (3S'(T,p ), Q(p) = a<p(T, p) — /3S'(T, p). Then 

D(p) = d(p) + (l + ap), (5) 

V ( r, p) = T (d(p) + Q(p)), S'(T, p) = ± (d(p) - Q(p)). 

Since <p(x, p)S'(x, p) — (p'(x, p)S(x, p) = 1, it follows that 

g 2 (p) = D 2 (p) - 4a/3(l + <p'(T, p)5(T, p)), (6) 

and consequently, 

Q(p)Q(p) = i>W£>W - 2a/3(<p’(T, p)S(T, p) + p’(T, p)S(T, p)), (7) 


where ’’dot” denotes derivatives with respect to p. Let n e I. Denote 



Q(Vn) = 0, 

QM 0, arg Qiy n ) e [0,7r), 
Q{v n ) 0, arg Q{vn) G [ir, 2n), 


0J r , 


:= d[ u \u n 


v = 0, m n — 1, Iq = {n G /' 




= 0}, I\ = {n G I' 


sequence D = {c o n } n &i U {Unu} n ei 0 v= \,m n -i i s called the D - sequence for L. 


r = 


i.e. m. 


= 1 for all n ), then D = {o ; n } ne a 


uj n ^ 0}. The 
We note that if 


Let a,/3 and q j are known a priori and fixed. The inverse problem is formulated as follows. 
Inverse problem 1. Given a(p),d(p) and D, construct B. 

Obviously, in general it is not possible to recover all coefficients from (2)-(3). Note that 
this inverse problem is a generalization of the classical inverse problems for Sturm-Liouville 
operators [4-7]. 


3. Properties of the spectral characteristics. 


Let $(a;, p) be the solution of equation (1) under the jump conditions (3) and the boundary 
conditions $(0,p) = 1, $(T,p) = 0. Denote M(p) := &(0,p). The function M(p) is called 
the Weyl-type function. Clearly, 


p) = C(x, p) + M(p)S(x, p), 


(8) 
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M(P) = -ff (9) 

Since {C(x, p), S(x, p)) = 1, it follows from (8) that 

{Hx,p),s(x,p)) = l. ( 10 ) 

Denote T k := b k — b k - 1 , k = 1,1V. Then = 7\ + • • • + T k , T = 7\ + ... + T N . Put 
x k = x - b k _ i for x G [b k -i,b k ]-, hence x k G [0 ,T k ], Let S k (x k ,p ) and C k (x k , p) be the 
solutions of equation (1) on [b k -i,b k ] under the initial conditions 

S k {0,p) = c' k ( 0 ,p ) = 0 , S k ( 0 , p) = C k (Q, p) = 1 . ( 11 ) 

For each hxed x k , the functions S k \x k , p) and C k '\x k , p), v = 0,1. are entire in p, and 

(C k (x k , p), S k (x k , p)) = 1. 

Lemma 1. The following relations hold for k = 1, N — 1, v — 0,1 : 

s w (fc+1 - o,p) = 7 kS(h - o ,p)c { ^(t m ,p) + 7 ps\w - o, P )s { th(n + up) 

+(*M + - 0, p)S^i(T t+ i, p), (12) 

C ( ->(6 1+1 - 0,p) = 7- O.pJC&pi+i.p) + T pC'(&* - O.pJS&pUi.p) 

+(*P'?fc + Vi)C(h - 0,p)52,(Ti + i,p), (13) 

Indeed, £x k = 1,iV — 1. Let x G [6fc,6fc+i], i.e. x = x k+ \ + b k , x k+ \ G [0, T k+1 ], Using 
the fundamental system of solutions C k+ i(x k+ i, p), S k+ i(x k+ i, p), one has 

S {u) (x,p) = A(p)Ci%(x k+1 , p) + B(p)S^ 1 (x k+1 ,p), u = 0,l. 

Taking initial conditions (11) into account we find the coefficients A(p) and B(p), and arrive 
at (12). Relation (13) is proved similarly. 

Denote 

1 f x 1 f T 

£( x ) = g J PW dt ’ u= fff I p (t) dt, t = Im p, G s = {p: \p - v n \ > 5 Vn}, 

n± = {p : ±t > 0}, Uf = {p : arg p G [S,n - 5]}, Uj = {p : arg p G [v: + 8, 2vr - 5]}. 

It is known (see [9]) that there exists a fundamental system of solutions 
{Y? (x,p), Y±(x,p)}i x G [0, T], p G Ll^, of equation (1) with the properties: 

1) The functions Yff(x, p) are regular in p G IP* 1 , |p| > p*, and are continuous for x G 
[0, T], p G II ± , |p| > p* . 

2) For \p\ —> 00 , p G II± k = 1, 2, z/ = 0,1, 

^ Y k( x ’p) = {p R k) u exp(px + £(x))R k )[l], [1] = l + 0(p _1 ), = (-l) fc_1 L (14) 

Using (14) and Lemma 1, one gets for x G (bj,bj+ 1 ), |p| —* 00 , p G II ^ : 

£± _ _ £± 

C (!/) (x,p) = 1 2 J (T^p) i/ exp(=Ff(p^ + g(^))[l], (15) 

S {v \x,p) = T ' 2 . 3 (Tip)" e*p(Ti(px + £(a;))[l], (16) 

® M (hP) = c ± 1 . ± (±^p) iy exp(±i(pa: + g(a:))[l], (17) 

£1 •••U 



4 


In particular, we have for \p\ —y oo, p G II^ : 

«(P) = y (f* ''' ^-i) ex P(T^(p+u;)T)[l], d(p) = T^(^ 1 ± ...^_ 1 )exp(Y(p+w)r)[l]- (18) 
Moreover, for x G [0, T], p G LU : 

|S H 0r,p)| < C\p\ v ~ l exp(|r|x), \C {u \x,p)\ < C\p\ u exp(\r\x), (19) 

I^Hup)! < C\p\ v exp(~\r\x), \M(p)\ < C\p \, p e G s . (20) 

Let n G I. Using (9) we obtain that in a neighborhood of the point p = v n , the function 
M(p) has the representation 


M(p) = 


m n — l 

^ Ip 

v =0 


M, 


n+v 


Vr, 


W+l 


+ M* (p), 


( 21 ) 


where M*(p ) is regular in p — u n , and the coefficients M n+V , v = 0, m„ — 1 are calculated 
from dl\is n ) and S u+mn ' ) (u n ) for u = 0, m n — l. More precisely, 


M, 


n+m n — 1 —v 


1 1 

~J i^un ”b ^ ^ I 'd n -\-m n — l—kdp—k,ri \ j ^ 0, TTl n 1, 

«0n ' ' 


fc =0 


p 1 •— \ p — _ 1 _ 

u i/7i I U 1 vnp “vre / . \| 

zd (z/ + m„)! 


d(* /+m » ) (i/ n ), z/ = 0, m n — 1. 


In particular, M n+m „_i = —df n /do n . If w-n = 1 (i.e. n E I \ I'), then 


di(r'n) 


M n = 

d(i'n) 


d(p) := -j- d(p). 

dp 


( 22 ) 


The sequence {M n } nS A is called the Weyl sequence. The data X> = {z/ n , M n } nS A are called 
the spectral data. We note that the specification of the spectral data X> uniquely determines 
the Weyl-type function M(p) (see [16]). 

We consider the following auxiliary inverse problem which is called IP-0. 

IP-0. Given the spectral data V = {u n ,M n } ne a and 7 j, j = 1,1V — 1, construct 
P(x),q(x),xe (0,T), p'jppj, j = l,N-l. 

Let us prove the uniqueness theorem for IP-0. For this purpose together with B we consider 
a boundary value problem B of the same form but with different potentials p(x),q(x), and 
different coefficients of the boundary and the jump conditions. We agree that if a sertain symbol 
6 denotes an object related to B , then 6 will denote the analogous object related to B. 

Theorem 1. If V = V, 7 j = 7 j, j = 1,1V — 1, then p = p,q = q, 77 ' = = fjj , 

J = 1,1V — 1. 

Proof. Consider the functions 


Pn(x,p) = $(x,p)S'(x,p) - S(x,p)&(x,p), 1 
Pn(x,p) = S(x,p)$(x,p) - ${x,p)S(x,p). J 


(23) 


Using (10), we calculate 


S(UP) = Pn(x,p)S(x,p) + Pi 2 (x,p)S'(x,p), 1 
<F(a:,p) = Pn(x, p)$(x, p) + Pi 2 (x, p)&(x, p). J 


(24) 
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According to (8) and (9), for each fixed x, the functions Pn(x,p) and Pi 2 (x,p) are mero- 
morphic in p with poles at the points v n and v n . Denote = G$ D G$■ By virtue of (19), 

(20) and (23) we get 


\Pn(x,p)\<C\p\-\ \P u (x,p)\<C, pe&l (25) 

It follows from (8) and (23) that 

Pn(x,p) = C(x,p)S'(x,p ) - S(x,p)C'(x,p) + (M(p) - M(p))S(x,p)S'(x,p), 

Pi 2 (x,p) = S(x, p)C(x, p) - C(x, p)S(x, p) - (M(p) - M(p))S(x, p)S(x, p). 

Since V = T>, it follows from (21) that for each fixed x, the functions Pik(x,p) are entire in 
p. Together with (25) this yields Pu(x,p) = 0, Pn(x,p) = A{x), where the function A{x) 
does not depend on p. Using (24) we derive 

S(x,p ) = A(x)S(x,p), Q(x,p) = A(x)$(x, p). (26) 


Taking (15)-(17) and (23) into account, we obtain for each fixed x G (bj,bj + 1 ) : 

+ vr—-r)[i], pe n^, |p| ->■ oo. 


Pll(x,p) = -y~ ± 


1^7.C; , .C; 






Therefore, the function A(x) is piecewise constant (a step-function). Together with (26) this 
yields q(x) = q(x),p(x) — p(x), x G (0,T). It follows from (26) that 


S{x,p) = ${x,p) 

S(x,p ) ${x,py 

and consequently, (£f) 2 = {£,f) 2 i j — 1,N — 1. Then r/' = f/' , j = 1, iV — 1, and con- 
seq uently, = §= , j = 1,77-1, S(x,p) = S(x,p ), &(x,p) = $(x,p), rjj = fjj , 
j — 1, N — 1. Theorem 1 is proved. 

Using the method of spectral mappings [3], one can obtain a constructive procedure for the 
solution of IP-0 (see [3] for details). 


4. Solition of the Inverse problem 1. 


Let a(p),d(p) and be given. Note that a, /3, 'jj, j — 1, N -£;1 are known a priori. 

The solution of the Inverse problem 1 is constructed as follows. 

First we calculate the zeros V = (z / n } ne a of d(p). Using the asymptotics (18) we find 
4 ) and h 1 . Taking (5) into account, we construct D(p). According to (6) we 
calculate Q 2 [y n ) = D 2 {y n ) —4a/3, and 

Q{v n ) = Un\/D 2 {v n ) - Aa(3, n G I. (27) 

Here and below we agree that if z — \z\e^, £ G [0,27r), then yT = |z| 1/,2 e^/ 2 . We construct 
ujno by 

&n0 — ^(-^(Ui) + Q( U rSji (28) 

since co n0 = di(u n ) = C(T, u n ) = p(T, v n ). 

We construct the Weyl sequence {M n } ne \ as follows: 

Case 1. Let n G / \ /' (i.e. m n = 1 ). Then, in view of (22), 


d{y n ) 


(29) 
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Case 2. Let n E h (i.e. m n >1, co n fi 0 ). Then it follows from (7) that 

(Q(p)Q(p))\pZll = (D(p)D(p))\ u p zll u = hrh~^i. (30) 

Using (30) we find Q^\u n ), v = l,m n — 1. Since 

di(p) = ^( D (p) + Q(P)) ~ (ipti + h)d(p), 


we construct u nv , v = l,m n — 1 by the formula 


u nv = ( D M (is n ) + (z/ n ) j, v = 1, m n - 1. (31) 

Case 5. Let n e J 0 (i.e. m n > 1, o; n = 0 ). Then cn ni/ , u = l,m n — 1 are given a priori. 

Thus, we have constructed the Weyl sequence {M n } ne By solving the auxiliary inverse 
problem IP-0 we find p(x),q(x),x G (0, T), r/', r/j, j = 1,N — 1, and then we calculate the 
coefficient /i, using (4). Thus, the following theorem holds. 

Theorem 2. T/ie specification of a(p),d(p) and uniquely determines the boundary 
value problem B. The solution of Inverse problem 1 can be found by the following algorithm. 
Algorithm. Given a(p),d(p) and 17. 

1) Calculate zeros V = {v n } n£ a of d(p). 

2) Find hi using (18). 

3) Construct D(p) by (5). 

f) Calculate Q 2 {v n ) = D 2 {y n ) —4 a/3, and Q(v n ) by (27). 

5) Construct oj n o by (28). 

6) Construct Q^(u n ), n Eli, v — 1 ,m n — 1 using (30). 

7) Find io nv , n Eh, v — 1 ,m n ~ 1, via (31). 

8) Calculate the Weyl sequence {M n } n >i using (29) and the recurrent formula 


M, 


n+rrin — 1—v 


l 

don 


u-1 


ifi^vn ^ hWn+m n — l—kdv—k,ri') , Tl E I, V 0, 'Tfl n 1, d un . 


' W. L 


k =0 




9) Find p{x),q{x),x E (0,T) and q'jPPj, j — 1, N — 1, by solving the inverse problem IP-0. 

10) Calculate the coefficient h . 

Similarly, one can solve the following inverse problem. 

Inverse problem 2. Given V, V and Tl, construct B. 
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